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Abstract — A validation study of Computational Fluid Dynamics 
(CFD) for Supersonic Retropropulsion (SRP) was conducted 
using three Navier-Stokes flow solvers (DPLR, FUN3D, and 
OVERFLOW). The study compared results from the CFD codes 
to each other and also to wind tunnel test data obtained in the 
NASA Ames Research Center 9'x7' Unitary Plan Wind Tunnel. 
Comparisons include surface pressure coefficient as well as 
unsteady plume effects, and cover a range of Mach numbers, 
levels of thrust, and angles of orientation. The comparisons 
show promising capability of CFD to simulate SRP, and best 
agreement with the tunnel data exists for the steadier cases of 
the 1-nozzle and high thrust 3-nozzle configurations. 
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1. Introduction 

Supersonic Retropropulsion (SRP) is a potentially viable 
means to decelerate high mass vehicles entering the Martian 
atmosphere [l]-[6]. Previous methods of deceleration are not 
scalable for exploration type vehicles which can potentially 
weigh tens of metric tons. Since ground and flight testing of 
SRP at entry conditions can be difficult and cost-prohibitive, 
the development of this enabling technology can be enhanced 
with the ability to predict the flowfield numerically using 
Computational Fluid Dynamics (CFD). 

SRP results in a complex flow structure involving shocks, 
shear layers, recirculation and stagnation regions, which 
makes validation of the CFD methods a high priority. The 
validation process includes using multiple CFD codes to com- 
pare to historical and recent wind tunnel tests. Three CFD 
codes are being applied to SRP: DPLR (Data Parallel Line 
Relaxation) [32], FUN3D (Fully Unstructured Navier-Stokes 
Three-Dimensional) [3 3] [34], and OVERFLOW (OVERset 


grid FLOW solver) [35]. The codes all solve the Navier- 
Stokes equations, but differ in implementation, grid type, 
and numerical methods. Through code-to-code and code-to- 
test comparisons, best practices in grid generation, numerical 
method selection, and solution advancement are established, 
and validity is added to the CFD methods. With continuing 
validation, confidence is built for using CFD to predict Mars 
entry conditions. 

CFD validation efforts for SRP were conducted under the 
NASA Exploration Technology Development Program. The 
CFD solvers were applied [8] to a historical [7] and a more 
recent [36] wind tunnel test. Although much was learned 
through the exercise, the historical test reports lacked key 
information for CFD validation. A wind tunnel test was then 
conducted for the purpose of CFD validation in the Langley 
Research Center supersonic 4' x 4' Unitary Plan Wind Tunnel 
(LRC UPWT) [14], the design of which was aided by CFD 
[12] [13]. Through results from the test, the CFD best prac- 
tices were established [15] and an extensive comparison study 
was conducted [16]. The study yielded promising results, 
but not all conditions were properly predicted consistently 
between the CFD codes. To further the validation exercise, 
and to obtain higher and more flight-relevant thrust levels, the 
same model was tested in the NASA Ames Research Center 
9' x 1' Unitary Plan Wind Tunnel (ARC UPWT). This paper 
will focus on the subsequent CFD validation study. A subset 
of the study can be found in reference [17], upon which this 
paper builds. 

Qualitative comparisons of the flow structure will be made 
by comparing CFD to high-speed shadowgraph images, and 
quantitative comparisons will be made by comparing aver- 
aged surface pressure with pressure tap data from the tunnel. 
High frequency pressure transducer data are also available, 
and have been studied more thoroughly for the LRC UPWT 
test [37] than for the ARC UPWT test. Comparisons of 
dominant frequencies of certain runs were discussed in ref- 
erences [15] and [17]. Neither test used a balance system, so 
forces and moments were not measured. However, code-to- 
code comparisons of forces and moments were included in 
references [16] and [17]. 

This paper will first introduce the SRP flow structure, the 
CFD solvers, and the ARC UPWT test. It will then dis- 
cuss the different modes the SRP flow structure experiences 
depending on nozzle configuration, thrust level, angle of 
attack, and Mach number. Then CFD results are compared 
to tunnel data, and strengths and weaknesses of the solvers 
are discussed. 
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2. SRP Flow Structure 

Figure 1 shows a CFD simulation of a run from the Daso, et 
al. wind tunnel test [36] to represent the SRP flow structure. 
In a supersonic freestream, a bow shock forms around the 
model, which is an Apollo capsule with an embedded sonic 
nozzle. As an opposing jet is initiated, the bow shock is 
pushed upstream as the apparent body frontal area is in- 
creased with the appearance of a barrel plume. The barrel 
plume contains free shear layers as well as a terminal shock. 
Between the terminal and bow shocks is an interface or 
contact surface where opposing streams stagnate. With the 
barrel plume, recirculation regions appear as well as a triple 
point, where three types of flow meet: supersonic jet flow, 
subsonic post-shock jet flow, and subsonic recirculating flow. 
The term “triple point” has a two-dimensional sense to it, in 
three dimensions it is an annular ring. 

This is an example of 1 -nozzle SRP flow at a relatively high 
thrust level. Other modes exist depending on the jet-to- 
freestream pressure ratios, including a long-jet penetration 
mode at low pressure ratios [7]. Even other modes exist 
and can vary according to thrust level, angles of attack, 
and multiple nozzle configurations where plume-to-plume 
interactions may exist. 


3. Description of CFD Solvers 

The three solvers applied to the SRP problem differ in 
implementation, grid type, and numerical methods. DPLR 
and FUN3D are finite- volume while OVERFLOW is finite- 
difference. DPLR uses cell-centered structured overset grids, 
while OVERFLOW uses node-centered structured overset 
grids. FUN3D employs node-centered unstructured grids. 
With these differences between codes, much is to be learned 
through code-to-code comparison when applying them to a 
single set of problems. 

In reference [15], grid and time resolution, observed order 
of accuracy, and the establishment of best practices are de- 
scribed in detail. The best practices were established through 
CFD sensitivity studies being applied to a 1 -nozzle case from 
the LRC UPWT test. The following are brief descriptions of 
each code including the numerical methods used by each in 
this study. 

DPLR 

The DPLR CFD code is a parallel, structured multi-block, 
finite-volume code with overset grid capability that solves 
the Reynolds- Averaged Navier- Stokes (RANS) equations for 
continuum flow, including finite-rate chemistry and thermal 
non-equilibrium. In the present study, the thermally- and 
calorically-perfect RANS equations for air are solved implic- 
itly with lst-order time accuracy. Inviscid fluxes are formed 
via a modified Steger- Warming flux vector splitting [21] with 
3rd-order Monotone Upwind Scheme for Conservation Laws 
(MUSCL) extrapolation [22] subject to a minmod limiter 
[23] and 2nd-order flux integration. The viscous terms are 
computed with a 2nd-order spatial accuracy central difference 
approach. For the present analysis, Menter’s Shear-Stress 
Transport (SST) turbulence model [24] was employed with 
a vorticity-based production term and no compressibility 
corrections. 

FUN3D 

FUN3D contains a node-based finite- volume unstructured 
flow solver. For this study, FUN3D was run with a 
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Figure 1 . 1 -nozzle SRP flow structure description diagrams. 


selectively-dissipative version of the Low-Dissipation Flux 
Splitting Scheme (LDFSS) inviscid flux function [25] and 
a modified Van Albada [26] limiter according to Vatsa and 
White [27]. Direct Eddy Simulation (DES) turbulence mod- 
eling was used with Spalart-Allmaras (SA) [28] as the sub- 
model. Node-based conservative variables are computed by 
driving a 2nd-order accurate spatial residual to steady-state 
with a point-implicit iterative method. A modified, optimum 
2nd-order backward difference formula (BDF) scheme is 
used in conjunction with a temporal error controller that 
assured design order [29]. 
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OVERFLOW 


OVERFLOW 2 is an implicit Navier-Stokes flow solver that 
uses structured overset grids. For the current work, the 
HLLE++ numerical flux function [30] with the Van Albada 
limiter was used for spatial terms, and the Symmetric Suc- 
cessive Over Relaxation (SSOR) algorithm with dual time 
stepping using Newton sub-iterations for temporal terms. 
DES turbulence modeling was used with SST as the sub- 
model, the SST model used the strain-based production term 
and employed Wilcox’s realizability constraint [31]. The 
overall scheme is 2nd-order accurate in space and time. The 
calculation of the inviscid fluxes for both the flow solver 
and the turbulence model use 3rd-order accurate MUSCL 
reconstruction and 2nd-order flux quadrature. 


4. SRP Wind Tunnel Test 

Test 234 in the ARC UPWT was designed specifically for 
SRP CFD validation and was a follow-on to test 1853 of the 
LRC UPWT. The reason for the ARC UPWT test was to use 
the same model from the previous test in a larger test section 
to be able to obtain higher and more flight-relevant thrust 
coefficients (Cr = T/(qA)). 

The model was a 70° sphere-cone forebody with a 5 inch 
diameter cylindrical side body. The model included four 
nozzles which could be plugged to offer 0-, 1-, 3-, or 4-nozzle 
configurations. One nozzle was located at the center of the 
forebody, and the three others were oriented radially every 
120° at the model’s 1/2 radius. Air was used as both the 
freestream and jet gases. The test data included high speed 
shadowgraph movies (5-10 thousand frames per second) and 
pressure readings from 167 taps including seven 40 kHz pres- 
sure transducers. From these data, qualitative comparisons 
with CFD flow structure and unsteady behavior can be made 
with the shadowgraph movies, as well as averaged surface 
pressure comparisons with the pressure taps. 

The test run matrix included four nozzle configurations: 0, 1 
(center), 3 (peripheral), and 4 (center and peripheral). Two 
Mach numbers were tested: 1.8 and 2.4. The Reynolds 
number per foot for Mach 1.8 was 1.5 million, and for Mach 
2.4, 1 and 1.5 million were tested. The total freestream 
temperature for Mach 1.8 was 550 degrees Rankine, and for 
Mach 2.4 was 555 degrees Rankine. For angles of orientation, 
a ranged from -8° to 20° and /3 ranged from 0° to 12°. 
Thrust coefficients ranged from 0 to 10 in increments of 2. 
To mitigate liquefaction, the jet air supply was heated to 
150° F; however, evidence of some liquefaction was noted 
at high thrust levels. A full description of the test can be 
found in reference [20]. Figure 2 is a snapshot of the model 
in the 1 -nozzle configuration in the test section. Figure 3 
shows diagrams of the model face and side-body pressure tap 
layouts [19]. 

Tunnel uncertainty will not be included in this paper since 
calculations are still underway, but will be carried out with 
the same method as the LRC UPWT test [19] as described in 
reference [9]. 



Figure 2. SRP test model in the 1 -nozzle configuration in 
the NASA ARC 9' x 7' UPWT. 


Forebodv Instrumentation 



o ESP Ports 
• Kulites 


Aftbodv Instrumentation 
x/L 


0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 



Figure 3. Pressure tap layout on the test model. Top image 
is the front face of the model and the bottom image is the 
side body. Filled circles represent the high frequency Kulite® 
pressure transducers. Image represents 172 planned pressure 
taps including nine 40 kHz pressure transducers. The final 
count was 167 taps with seven 40 kHz pressure transducers. 
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5. SRP Flow Structure Modes 

The purpose of this paper is to explore the different modes 
that exist for SRP when varying Mach number, the number of 
jets, the levels of thrust, and the angle of orientation, and how 
well CFD can match those modes in behavior and in average 
surface pressure. SRP literature shows that different modes 
exist for the 1 -nozzle configuration depending on the jet-to- 
freestream pressure ratio [7], [36]. At low ratios, a long-jet 
penetration mode exists where the jet exhaust is unsteady, 
similar to a candle flame, and the bow shock has a large and 
fluctuating standoff distance from the model. At higher ratios, 
the short-penetration mode exists (described in Section 2) 
where the flow is much steadier and the bow shock has a 
smaller standoff distance than the long-jet penetration mode. 
It was observed in the LRC UPWT test that once in the short- 
penetration mode, unsteadiness in the triple point can vary 
between a periodic oscillation and a smaller scale aperiodic 
fluctuation. Also, as angle of attack increases, the windward 
triple point oscillations become larger as the leeward triple 
point oscillations decrease, then for a = 20° the steadiness 
of the shock structure was lost as the barrel plume and bow 
shock fluctuated chaotically. The CFD solvers did well 
simulating the 1 -nozzle cases of the LRC UPWT test in 
behavior and in matching the average surface pressure, with a 
trend that DPLR predicted a steadier flowfield than was seen 
in the test, and FUN3D and OVERFLOW well predicted the 
unsteady behavior seen in the test. The different levels of 
predicted unsteadiness in the codes is attributed to turbulence 
modeling and grid resolution [15], [16]. 

More observations from the LRC UPWT test showed that 
for the 3 -nozzle configuration at lower thrust coefficients 
(< 3), the jet-to-jet interactions and bow shock behavior 
were aperiodically unsteady. The CFD did well in qualitative 
comparisons of the shock structure behavior, and varying 
success was obtained when comparing to the average surface 
pressure. Discrepancies between the codes and test data at 
the model nose and shoulder implied a difference in how 
the codes were simulating the jet-to-jet interaction and the 
unsteady bow shock shedding and how those impacted the 
model face. A large discrepancy in data acquisition rates did 
exist between the CFD and wind tunnel test (WTT). Rates 
were on the order of 10 Hz for the WTT and 10 kHz for the 
CFD simulations. An entire CFD run could fit within the time 
between two pressure tap readings, meaning high frequency 
fluctuations captured by the CFD codes were not captured 
by the WTT system, and the low frequency fluctuations 
captured by the WTT system were not captured by the CFD. 
This introduced an error source that increased when unsteady 
frequencies became very high or very low. A few solitary 
runs at zero angles of attack and Cp up to six showed that 
the bow shock became nearly steady as jet-to-jet interactions 
decreased due to distinctly defined barrel plumes. The more 
steady flowfield was desirable for vehicle stability, and this 
helped build the case for the ARC UPWT test where higher 
thrust was possible. As explained in reference [17], at higher 
thrust for the 3 -nozzle configuration (> 6) at zero angles of 
attack the flow was much steadier and the CFD compared 
better qualitatively and quantitatively than for the lower thrust 
cases from the LRC UPWT test. 

For the 4-nozzle configuration, the LRC UPWT test showed a 
short-penetration bow shock behavior at zero angle of attack. 
A discrepancy between the CFD codes for this behavior 
was noted. Runs 307 and 311 were the same in Mach 
number (4.6) and thrust coefficient (2), but differed in a 
roll angle of 0° and 180°, respectively. The test showed 
the same short-penetration shock behavior for both runs, 


but DPLR showed a larger standoff distance for run 307, 
and OVERFLOW showed a larger and fluctuating standoff 
distance for run 311. This implied a difference in modes 
for the 4-nozzle configuration, and the boundary between 
the modes was not well defined by the CFD codes. As 
discussed in [17], at higher thrust coefficients, the 4-nozzle 
configuration did demonstrate a large and fluctuating bow 
shock standoff distance, and became chaotically unsteady. In 
the LRC UPWT test, it was shown that for a roll angle of 
0°, unsteadiness increased with angle of attack. For a roll 
angle of 180°, angle of attack 12° was unsteady while angle 
of attack 20° settled into a steadier mode. All these behaviors 
were simulated well with the CFD solvers with the exception 
at angle of attack 0° already discussed. The average surface 
pressures for the LRC UPWT test compared well between the 
codes for all cases, with the largest discrepancies found at the 
model shoulder. Similar to the 3-nozzle configuration, the 
differences in average Cp at the shoulder implies differences 
in simulated bow shock shedding that impacts the model face. 


6. Results 

The CFD run matrix (Table 1) consisted of 3 1 cases varying in 
Mach number, nozzle configuration, thrust coefficient (Cp), 
and angles of orientation. All runs in the CFD matrix 
had a Reynolds number per foot of 1.5 million with the 
exception of run 223, which had 1.0 million. For the angles 
of orientation, most cases vary only in (3 angle since the [3 
plane provided better shadowgraph imagery for qualitative 
comparisons. Runs 279 and 309 included non-zero a angles 
and were meant to compare total alpha angles (ap) of 12° 
at different roll angles ( 0 ). Each nozzle configuration will 
now be discussed in more detail. The 0-nozzle configuration 
will be skipped since it was already covered in reference [17]. 
Also discussed in reference [17] are the effects of the sting on 
pressure measurements from the test. The sting caused an 
increase in the average Cp over the aft 30% of the model. 
None of the CFD simulations in this paper include the sting, 
so will underpredict the average Cp in that region. 

All cases in the run matrix were computed with OVERFLOW, 
while a smaller subset was computed with FUN3D and 
DPLR. For the flowfield images, unless otherwise labeled, 
the CFD solver is OVERFLOW. The line plots in this section 
compare the average surface pressure from the CFD results 
and the WTT data along 0 = 0° and 180° (see Figure 3). 
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Table 1 . CFD Run Matrix 


Run 

Mach 

Jets 

Cp 

a 

0 

4> 

106 

1.8 

1 

4 

0° 

0° 

180° 

116 

2.4 

1 

4 

0° 

0° 

180° 

223 

2.4 

1 

4 

0° 

0° 

0° 

223 

2.4 

1 

4 

0° 

-4° 

0° 

223 

2.4 

1 

4 

0° 

-8° 

0° 

223 

2.4 

1 

4 

0° 

-12° 

0° 

101 

1.8 

1 

2 

0° 

0° 

180° 

104 

1.8 

1 

3 

0° 

0° 

180° 

112 

2.4 

1 

2 

0° 

0° 

180° 

130 

1.8 

3 

6 

0° 

0° 

180° 

132 

1.8 

3 

8 

0° 

0° 

180° 

139 

2.4 

3 

4 

0° 

0° 

180° 

139 

2.4 

3 

4 

0° 

12° 

180° 

141 

2.4 

3 

6 

0° 

0° 

180° 

141 

2.4 

3 

6 

0° 

4° 

180° 

141 

2.4 

3 

6 

0° 

8° 

180° 

141 

2.4 

3 

6 

0° 

12° 

180° 

143 

2.4 

3 

8 

0° 

0° 

180° 

126 

1.8 

3 

2 

0° 

0° 

180° 

145 

2.4 

3 

10 

0° 

0° 

180° 

279 

2.4 

3 

6 

6° 

-10.33° 

30° 

309 

2.4 

3 

6 

-6° 

10.33° 

210° 

166 

1.8 

4 

2 

0° 

0° 

0° 

170 

1.8 

4 

6 

0° 

0° 

0° 

172 

1.8 

4 

8 

0° 

0° 

0° 

179 

2.4 

4 

4 

0° 

0° 

0° 

179 

2.4 

4 

4 

0° 

-4° 

0° 

179 

2.4 

4 

4 

0° 

-8° 

0° 

181 

2.4 

4 

6 

0° 

0° 

0° 

183 

2.4 

4 

8 

0° 

0° 

0° 

185 

2.4 

4 

10 

0° 

0° 

0° 


Figure 6 combines the results from all codes and the test for 
run 223, /3 = 0°. The DPLR simulation reached a steady 
state where FUN3D predicted a periodic oscillation of the 
triple point, and OVERFLOW predicted the more random 
triple point oscillation as seen in the test. On the model face, 
all codes overpredicted the average Cp near the shoulder, 
with FUN3D overpredicting the most due to simulating an 
organized shedding from the triple point. Along the model 
side, no significant differences were noted between the codes, 
and all underpredicted towards the aft end of the model due 
to not including the sting in the simulations. 


1 -Nozzle Configuration 

The highest Cp for the 1 -nozzle configuration in the ARC 
UPWT test was 4. The limiting factor was the high pressure 
system used for the jet exhaust. The CFD matrix focused on 
Cp of 4 with runs 106, 116, and 223, and how the behavior 
at that Cp varied with Mach number and (3 angle. Runs 
101, 104, and 112 looked into lower thrust coefficients. All 
cases demonstrated the short-penetration mode described in 
Section 2. 

Figure 4 visualizes the effect of Cp as seen in the test and 
how well the CFD matched that behavior. The increased size 
of the barrel plume was simulated well with OVERFLOW, 
and average surface pressures matched well, but consistently 
underpredicted in the aft region of the model due to not 
modeling the sting. For run 104, the average Cp was 
overpredicted on the model face, caused by OVERFLOW 
predicting a periodic oscillation of the triple point. Shedding 
from the triple point propagates to the model face, which in 
turn affects the Cp on the surface. When the shedding is 
periodic instead of random, the energy is more organized and 
the average Cp increases. Periodic shedding was not seen in 
the test for this run. 

The effects of increasing /3 angle is shown in Figure 5. As /3 
angle increased, the windward triple point oscillation became 
more pronounced as the leeward triple point became more 
stationary. This was predicted well by OVERFLOW, and the 
average surface Cp compared well on the model face with 
some underprediction on the model side. Run 223 /3 = 12° 
from the ARC UPWT test was similar in triple point behavior 
as run 165 a = 12° from the LRC UPWT test [15], [16]. 
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Figure 4. Effects of CV; test shadowgraph compared to CFD constructed Schlieren and average surface pressure comparison 
for the 1 -nozzle configuration at Mach 1.8. 
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Figure 5. Effects of /?; test shadowgraph compared to CFD constructed Schlieren and average surface pressure comparison 
for the 1 -nozzle configuration, run 223: Mach 2.4, Ct 4. 
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Figure 6. Code-to-code and code-to-test comparison of Run 
223: 1 -nozzle, Mach 2.4, Cp 4, (3 0°. 

Given the close comparisons in behavior and in average sur- 
face Cp, CFD fairs well for simulating the 1 -nozzle configu- 
ration. Some differences existed in the level of unsteadiness 
predicted by the solvers, which did not significantly affect 
differences in the average Cp. Two traits of the CFD that did 
affect the comparison to the tunnel data were sting effects and 
incorrectly predicting the triple point behavior. 

3-Nozzle Configuration 

As discussed in section 5, the higher thrust coefficients tested 
in the ARC UPWT demonstrated steadier flow than the lower 
thrust coefficients seen in the LRC UPWT test. This same 
trend is also shown in Figure 7 where run 126 was more 
unsteady at a Cp of 2 than runs 130 and 132 which were 
at Cps of 6 and 8, respectively. The OVERFLOW results for 
run 126 differed in behavior and in average surface Cp from 
the test, especially at the model nose. As Cp increased, the 
steadiness did as well, as did the agreement with the CFD 
simulations. Good agreement was seen between the codes, 
with the largest differences at the model nose and near the 
nozzle. The CFD underpredicted the Cp for the majority 
of the model face and side, with overpredictions at the nose 
by FUN3D for run 130 and by OVERFLOW for run 132. 
OVERFLOW closely matched the test average Cp at the nose 
for run 130, and near the nozzle for runs 130 and 132. DPLR, 
which reached a steady state solution, closely matched the 
average Cp on the model side for run 130. 

As (3 increased (see Figure 8) the test showed the shock struc- 
ture maintaining its basic shape, with a noticeable difference 
in the windward triple point of the windward plume which 


increased in oscillation size. This same behavior existed for 
all tested (3 angles up to 12°. Here lies a weakness in the 
OVERFLOW results. For (3 angles up to 8°, the windward 
plume maintained its shape, but for f3 = 12°, the windward 
plume was compressed, and the resulting flow structure was 
more unsteady than what was seen in the test. The difference 
in flow behavior did not greatly influence the average surface 
pressure comparison, all cases still compared well to the test 
data. 

The effects of roll (fi) at high angle of attack is shown in Fig- 
ure 9. In the test, the three barrel plumes still maintained their 
shapes regardless of roll, but in the OVERFLOW results, only 
<p = 210° showed that behavior where the rest showed the 
compressed windward plume. For the cases with = 180° 
and 30°, the compressed plume was more directly in line with 
the oncoming tunnel flow, or more relatively windward than 
for the = 210° case. For the = 210° case, the two 
windward plumes were equally windward, and neither were 
as directly in-lined with the flow as the other cases. 

The high pressure at the nose of the model for the 3-nozzle 
configurations should not be taken to be resultant of on- 
coming freestream flow penetrating the space between the 
plumes. Instead, the high pressure region is caused by a 
recirculation zone created by entrapment of jet exhaust and 
entrainment from surrounding flow. Differences in surface 
Cp between the codes can be attributed to differences in 
modeling the recirculation caused by that entrainment and 
entrapment. 
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Figure 7. Effects of Ct\ test shadowgraph compared to CFD constructed Schlieren/shadowgraph and average surface 
pressure comparison for the 3-nozzle configuration at Mach 1.8. The CFD solver visualizations for runs 126, 130, and 132 
are OVERFFOW, DPER, and FUN3D, respectively. 
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Figure 8. Effects of /?; test shadowgraph compared to CFD constructed Schlieren/shadowgraph and average surface pressure 
comparison for the 3-nozzle configuration, run 141: Mach 2.4, Ct 6. The CFD solver visualization for /? = 0° is FUN3D, and 
for the other /3 angles is OVERFFOW. 
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Figure 9. Effects of </>; test shadowgraph compared to CFD constructed Schlieren and average surface pressure comparison for 
the 3-nozzle configuration at Mach 2.4, Ct 6. 
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Run 139 was used for a code-to-code comparison and is vi- 
sualized in Figure 10. In the test, the flow was more unsteady 
than shown in the CFD solvers. FUN3D was closest to the 
unsteady behavior, but eventually reached a steady state. This 
demonstrated that the boundaries between different modes of 
behavior for SRP, in this case the boundary between unsteady 
and much steadier flow as Cp increases, was not captured 
well by the CFD codes. For the average Cp, FUN3D 
overpredicted at the nose and bounded the test data on the 
model side, DPLR closely matched on the face and side with 
an underprediction near the nozzle, and OVERFLOW closely 
matched on the face but predicted nearly symmetric Cp on 
the side. 

Confidence in the ability of CFD to simulate 3 -nozzle SRP 
configurations is stronger for higher Cp s than for lower 
Ct s. The boundary between the unsteady lower Cps and 
the much steadier higher Cp s is not well predicted by CFD, 
but once the steadier mode is obtained, better agreement 
between the CFD solvers and the test data is achieved. The 
LRC UPWT CFD validation exercise showed that all solvers 
closely matched the wind tunnel test qualitatively in behavior, 
but differences in average surface Cp did exist. A possible 
error source for the CFD may be the grid system, which was 
optimized for 1 -nozzle SRP flow, but may not be optimized 



Figure 10. Code-to-code and code-to-test comparison of 
Run 139: 3-nozzle, Mach 2.4, Cp 4. 


for multi-nozzle flows with interacting shear layers. 

4-Nozzle Configuration 

Section 5 discussed that differences in the shock standoff 
distance as well as the level of unsteadiness changed based 
on thrust coefficient for zero angle of attack 4-nozzle cases. 
At lower Cps the flow was steadier with a short bow shock 
standoff distance, and at higher Cp s the flow became chaotic 
with large and oscillating shock standoff distances. This trend 
was noted in the ARC UPWT test results, which ranged from 
a Cp of 2 for run 166 to a Cp of 10 for run 185 of the 
CFD matrix. The CFD matched the behavior of the zero /3 
cases well, but not much value was placed in the average 
surface pressure comparisons. Since the flow was chaotic and 
the data acquisition rates varied so greatly between the CFD 
and WTT data, the averages obtained by the CFD and WTT 
differed greatly. The WTT data acquisition rates were also 
not sufficient to gather enough samples for this chaotic flow. 

As /? angle increased, the level of unsteadiness decreased in 
the test, as seen with run 179 in Figure 11. This behavior 
was not matched by OVERFLOW, which continued to predict 
a highly chaotic flowfield at f3 angles of 4° and 8°. This 
may be a limitation in the grid system, which was optimized 
for a 1 -nozzle configuration and then applied to multi-nozzle 
configurations. 
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Figure 11. Effects of /?; test shadowgraph compared to 
CFD constructed Schlieren/shadowgraph and average surface 
pressure comparison for the 4-nozzle configuration, run 179: 
Mach 2.4, Cp 4. The CFD solver visualization for /3 = 0° is 
FUN3D, and for the other /? angles is OVERFLOW. 


Code-to-code comparisons can be seen in Figure 12 for run 
166 (Mach 1.8, Ct 4). The test still demonstrated the stead- 
ier, short shock standoff distance of lower thrust coefficients. 
This behavior was matched by FUN3D and OVERFLOW, 
but a large standoff distance was predicted by DPLR. The 
average surface Cp by DPLR overpredicted FUN3D and 
OVERFLOW, both of which were closer to the test data. 
The FUN3D Cp data showed a more dynamic behavior than 
OVERFLOW and DPLR, a result of insufficient averaging. 

The confidence in the ability of the CFD solvers to simulate 
4-nozzle flow is not as strong as it is for the 1 -nozzle configu- 
ration. The change in the short-penetration standoff distance 
to the wildly chaotic behavior was predicted to be at different 
points between the CFD solvers. Due to the highly unsteady 
behavior of the 4-nozzle configuration at Cps greater than 
four, it was difficult to obtain reasonable average surface 
pressure comparisons between the CFD and the WTT. Also, 
steadier modes at angle of attack were seen for the 4-nozzle 
configuration which were not simulated with the CFD. 

Thrust Dominance 

For the SRP cases discussed, the total axial force (C a, total = 
Cp + C a, aero) was dominated by the thrust coefficient and 
the contribution of the aerodynamic axial force was very 



Figure 12. Code-to-code and code-to-test comparison of 
Run 166: 4-nozzle, Mach 1.8, Cp 4. 


small. This can be seen in the bar chart in Figure 13 where 
the gray section of the bars represent Cp and the blue sec- 
tions represent C a, aero- This raises the question to whether 
aerodynamic effects are negligible for these types of SRP 
flows. More research needs to be conducted to answer this 
question. Entry trajectories and angles as well as spacecraft 
scale need to be taken into consideration. A high frequency 
oscillation at the WTT scale will be a much lower frequency 
at the flight scale, and may have more contribution to vehicle 
stability. Normal force and the vehicle moments also need to 
be taken into consideration, especially for non-zero vehicle 
entry angles. 



Figure 13. Comparison of the contributions of thrust and 
aerodynamics to the total axial force for the ARC UPWT 
CFD cases run by OVERFLOW. 
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Computational Cost 

Table 2 shows a breakdown of the computational cost for each 
of the codes for a typical SRP run. This is not a perfect 
comparison since the numbers represented are either just 
from a single run or averages of multiple runs. All cases were 
computed on the Pleiades NASA Advanced Supercomputing 
cluster, an SGI Ice cluster, on either Nehalem or Westmere 
nodes. Since each code used different amounts of time, 
number of iterations, number of sub-iterations, and number 
of grid points, a common metric of CPU seconds per iteration 
per grid point between the codes is reported for comparison. 

Part of the requirement for these validation runs was to test 
the ability of the solvers to predict the SRP flowfields. For 
this reason, time-accurate runs were necessary because it 
is impossible to know beforehand if the unsteady effects 
played a large role in the aerodynamic effects. Costs for 
validation cases were large; however, for production runs of 
parametric studies or database generation, the computational 
costs should be able to drop significantly as either steady state 
runs or higher efficient numerical and time step processes are 
adopted. 


Table 2. Computational Cost Per Case 


Solver 

CPU 

Hours 

Iterations 

Grid 

Points 

CPU 

s/iter/gp 

FUN3D 

28,000 

39,500 

4.2e7 

6.1e-05 

DPLR 

44,500 

106,000 

5.3e7 

2.9e-05 

OVERFLOW 

35,039 

73,500 

8.5e7 

2.0e-05 


7. Conclusions 

A second wind tunnel test for validation of CFD in SRP 
conditions was conducted in the ARC UPWT 9' x 7' tunnel. 
The ARC UPWT test had a larger test section than the 
previous test, which took place in the LRC UPWT 4' x 4' 
tunnel. With the larger test section, higher thrust coefficients 
were obtained that closer match flight requirements. 

For the 1 -nozzle configuration, not much difference was 
observed between the higher thrust coefficients of 4 obtained 
in the ARC UPWT and the lower thrust coefficients of 2 
obtained in the LRC UPWT. The short-penetration shock 
interaction structure was observed, and oscillations at the 
triple point increased for the windward triple point as angle 
of attack was increased. The CFD solvers compared well 
to the 1 -nozzle cases in the CFD matrix, only varying in 
predicting periodic oscillations of the triple point when ran- 
dom oscillations were observed in the test. Average surface 
pressure comparisons were promising, building confidence in 
the ability of CFD to predict 1 -nozzle configuration SRP. 

For the 3-nozzle configuration, the unsteady behavior noted 
at lower thrust coefficients of 2 did not exist for higher thrust 
coefficients where the flow structure becomes more steady, 
especially in the bow shock. This behavior was predicted 
well with the CFD, and average surface pressures agreed well 
with the WTT. The area of greatest deviation between the 
codes was along the model shoulder and at the model nose. 
The differences at the nose are attributed to differences in 
entrainment, entrapment, and recirculation being simulated 
between the codes at the model nose, as well as any shear 
layer interaction that may be occurring between the plumes. 
At (3 and olt 12°, OVERFLOW did not predict the same be- 


havior as the WTT, predicting a more compressed windward 
barrel plume for most cases. This discrepancy in qualitative 
comparison did not seem to greatly influence the average 
surface pressure comparison. 

For the 4-nozzle configuration, a steadier mode with a short 
shock standoff distance existed for lower Ct s. As Ct was 
increased, the flow becomes highly unsteady with the bow 
shock standoff distance aperiodically oscillating between 
small and large standoff distances. The unsteady flowfield 
was predicted well with the CFD codes, but the average sur- 
face pressures were not comparable to WTT data, presumably 
due to large differences in data acquisition rates between the 
CFD and WTT and the large difference in time- averaging 
windows. At [3 4° and 8°, a steadier mode was noted 
which was not predicted by OVERFLOW, which continued 
to predict a largely chaotic flowfield. 

Given the comparisons in this validation study, confidence 
was built in the ability of CFD to simulate 1 -nozzle con- 
figurations at all /3 angles tested, and high thrust 3 -nozzle 
configurations at (3 < 8°. The ability to simulate the 
behavior seen in the high thrust 4-nozzle cases was also 
demonstrated, but average surface pressure comparisons were 
poor. The ability of the CFD codes to predict the exact 
boundary between steady and unsteady modes for 1-, 3-, and 
4-nozzle configurations was not demonstrated, and in fact 
offered contradictory results in both the LRC UPWT and 
ARC UPWT tests for the 4-nozzle a/ (3 zero configuration. 
The average surface pressure comparisons were the best for 
the 1 -nozzle configuration and for the steadier runs of the 
3- and 4-nozzle configurations. This correlates to steadier 
flowfields where the large difference in data acquisition rates 
between the CFD WTT did not greatly influence the averages 
obtained. Another possible error source may be the grid 
systems used for each code which were optimized for 1- 
nozzle SRP flow, not for multi-nozzle flows with shear layer 
interactions. 

Thrust dominance and computational costs were also men- 
tioned. Computational costs were high for validation but 
could be much less for production. 

Future work for CFD validation would include live rocket 
tests including startup transients, more flight-relevant calcu- 
lations including Mars atmospheric and rocket nozzle gases, 
and an application to flight- scaled spacecraft to better under- 
stand the influence of aerodynamics on stability and control. 
Future wind tunnel tests could include pressure measurement 
instrumentation along the model shoulder to help validate 
CFD in that area. More CFD research may be needed to 
properly define the boundaries between different behavioral 
modes, as well as an optimization of a grid system for multi- 
nozzle SRP flowfields. 


Acknowledgments 

Pieter G. Buning of NASA Langley Research Center, Hamp- 
ton, Virginia; Phillip C. Stuart, Thomas M. Booth, and 
Darby J. Vicker of NASA Johnson Space Center, Houston, 
Texas provided valuable guidance on gridding and solver 
best practices for OVERFLOW. William T. Jones of NASA 
Langley Research Center, Hampton, Virginia provided the 
GridEx/batchEx unstructured grid generation framework and 
helped with its application to generate FUN3D grids. Todd 
R. White and Andrew J. Hyatt of ERC inc., Moffett Field, 
California helped develop DPLR best practices and pro- 


12 



vided overset gridding assistance for DPLR cases. David 
A. Saunders of ERC Inc., Moffett Field, California devel- 
oped the code used for OVERFLOW and DPLR simulated 
Schlieren/shadowgraph visualizations. 

The authors would like to acknowledge the support of the 
Exploration Technology Development and Demonstration 
(ETDD) Program, which was managed at NASA-Glenn Re- 
search Center. The work documented herein was performed 
as part of ETDD’s Entry, Descent, and Landing (EDL) Tech- 
nology Development Project, which was managed at NASA- 
Langley Research Center and supported by NASA-Ames 
Research Center, NASA-Johnson Space Center, and the Jet 
Propulsion Laboratory. 


References 

[1] Braun, R. D., and Manning, R. M., “Mars Exploration 
Entry, Descent, and Landing Challenges,” Journal of 
Spacecraft and Rockets, Vol. 44, No. 2, 2007, pp. 310- 
323. 

[2] Steinfeldt, B. A., Theisinger, J. E., Korzun, A. M., Clark, 

I. G., Grant, M. J., and Braun, R. D., “High Mass Mars 
Entry, Descent, and Landing Architecture Assessment,” 
AIAA Paper 2009-6684, September 2009. 

[3] Zang, T.A., and Munk, M. M., et al., “Entry, Descent 
and Landing Systems Analysis Study: Phase 1 Report,” 
NASA TM 2010-216720, July 2010. 

[4] Edquist, K. T., Dyakonov, A. A., Korzun, A. M., Shidner, 

J. D., Studak, J. W., Tigges, M. A., Kipp, D. M., Prakash, 
R., Trumble, K. A., Dupzyk, I. C., “Development of 
Supersonic Retro-Propulsion for Future Mars Entry, De- 
scent, and Landing Systems,” AIAA Paper 2010-5046, 
June 2010. 

[5] Korzun, A. M., and Braun, R. D., “Performance Charac- 
terization of Supersonic Retropropulsion for High-Mass 
Mars Entry Systems,” Journal of Spacecraft and Rockets, 
Vol. 47, No. 5, 2010, pp. 836-848. 

[6] Korzun, A. M., Braun, R. D., and Cruz, J. R., “Survey of 
Supersonic Retropropulsion Technology for Mars Entry, 
Descent, and Landing,” Journal of Spacecraft and Rock- 
ets, Vol. 46, No. 5, 2009, pp. 929-937. 

[7] Jarvinen, P. O., and Adams, R. H., “The Aerodynamic 
Characteristics of Large Angled Cones with Retrorock- 
ets,” NASA CR NAS 7-576, February 1970. 

[8] Trumble, K. A., Schauerhamer, D. G., Kleb, W. L., Carl- 
son, J. R., Buning, P. G., Edquist, K. T., and Barnhardt, 
M. D., “An Initial Assessment of Navier-Stokes Codes 
Applied to Supersonic Retro-Propulsion,” AIAA Paper 
2010-5047, June 2010. 

[9] Aeschliman, D. P. and Oberkampf, W. L., “Experimental 
Methodology for Computational Fluid Dynamics Code 
Validation,” AIAA Journal, Vol. 36, No. 5, 1998, pp. 733- 
741. 

[10] Oberkampf, W. L. and T. G. Trucano, “Verification and 
Validation in Computational Fluid Dynamics.” Progress 
in Aerospace Sciences, Vol. 38, No. 3, 2002, pp. 209- 
272. 

[11] Oberkampf, W. L. and Roy, C. J., Verification and 
Validation in Scientific Computing, Cambridge, UK, 
Cambridge University Press, 2010. 

[12] Trumble, K. A., Schauerhamer, D. G., Kleb, W. L., 
and Edquist, K. T., “Analysis of Navier-Stokes Codes 


Applied to Supersonic Retro-Propulsion Wind Tunnel 
Test,” IEEEAC Paper 1471, March 2011. 

[13] Berry, S. A., Laws, C. T., Kleb, W. L., Rhode, M. N., 
Spells, C., Mccrea, A. C., Trumble, K. A., Schauerhamer, 
D. G., Oberkampf, W. L., “Supersonic Retro-Propulsion 
Experimental Design for Computational Fluid Dynamics 
Model Validation,” IEEEAC Paper 1499, March 2011. 

[14] Berry, S. A., Rhode, M. N., Edquist, K. T., and Player, 

C. J., “Supersonic Retropropulsion Experimental Results 
from the NASA Langley Unitary Plan Wind Tunnel,” 
AIAA Paper 2011-3489, June 2011. 

[15] Kleb, W. L., Carlson, J. R., Buning, P. G., Berry, 
S. A., Rhode, M. N., Edquist, K. T., Schauerhamer, 

D. G., Trumble, K. A., Sozer, E., “Toward Supersonic 
Retropropulsion CFD Validation,” AIAA Paper 2011- 
3490, June 2011. 

[16] Schauerhamer, D. G., Trumble, K. A., Kleb, W. L., 
Carlson, J. R., Edquist, K. T., “Continuing Validation 
of Computational Fluid Dynamics For Supersonic Retro- 
propulsion,” AIAA Paper 2012-0864, January 2012. 

[17] Zarchi, K. A., Schauerhamer, D. G., Kleb, W. L., 
Carslon, J. R., Edquist, K. T., “Computational Fluid Dy- 
namics Validation and Post-test Analysis of Supersonic 
Retropropulsion in the Ames 9x7 Unitary Tunnel,” AIAA 
Paper 2012-2705, June 2012. 

[18] Rhode, M. N., Oberkampf, W. L., “Estimation of 
Uncertainties for a Supersonic Retro-Propulsion Model 
Validation Experiment in a Wind Tunnel,” AIAA Paper 
2012-2707, June 2012. 

[19] Rhode, M. N., and Oberkampf, W. L., “Estimation of 
Uncertainties for Supersonic Retropropulsion Test 1853 
NASA Langley Research Center Unitary Plan Wind Tun- 
nel,” NASA TR-EDL_0 1 _TR_92 1 7 , October 2011. 

[20] Berry, S. A., Rhode, M. N., Edquist, K. T., “Supersonic 
Retropropulsion Experimental Results from the NASA 
Ames 9- x 7-Foot Supersonic Wind Tunnel,” AIAA Paper 
2012-2704, June 2012. 

[21] MacCormack, R. W. and Candler, G. V., “The Solution 
of the Navier-Stokes Equations Using Gauss-Seidel Line 
Relaxation,” Computers and Fluids, Vol. 17, No. 1, 1989, 
pp. 135-150. 

[22] van Leer, B., “Towards the Ultimate Conservative 
Scheme. V. A Second-Order Sequel to Godunovs 
Method,” Journal of Computational Physics, Vol. 32, 
1979, pp. 101-136. 

[23] Yee, H. C., “A Class of High-Resolution Explicit and 
Implicit Shock Capturing Methods,” NASA TM 101088, 
February 1989. 

[24] Menter, F. R., “Two-Equation Eddy- Viscosity Turbu- 
lence Models for Engineering Applications,” AIAA Jour- 
nal, Vol. 32, No. 8, 1994, pp. 1598-1605. 

[25] Edwards, J. R., “A Low-Diffusion Flux-Splitting 
Scheme for Navier-Stokes Calculation,” Computers and 
Fluids, Vol. 26, No. 6, 1997, pp. 653-659. 

[26] van Albada, G. D., van Leer, B., and Roberts, W. W., “A 
Comparative Study of Computational Methods in Cosmic 
Gas Dynamics,” Astronomy and Astrophysics, Vol. 108, 
No. 1, 1982, pp. 76. 

[27] Vatsa, V. N. and White, J. A., “Calibration of a Uni- 
fied Flux Limiter for Ares-Class Launch Vehicles from 
Subsonic to Supersonic Speeds,” JANNAF Paper, April 
2009. 


13 



[28] Spalart, P. R. and Allmaras, S. R., “A One-Equation 
Turbulence Model for Aerodynamic Flows,” AIAA Paper 
92-0439, January 1992. 

[29] Vatsa, V. N., Carpenter, M. H., and Lockard, D. R, 
“Re-evaluation of an Optimized Second Order Backward 
Difference (BDF20PT) Scheme for Unsteady Flow Ap- 
plications,” AIAA Paper 2010-0122, January 2010. 

[30] Tramel, R., Nichols, R., and Buning, R, “Addition 
of Improved Shock-Capturing Schemes to OVERFFOW 
2.1 ” AIAA Paper 2009-3988, June 2009. 

[31] Wilcox, D. C., “Formulation of the k - uo Turbu- 
lence Model Revisited,” AIAA Journal, Vol. 46, No. 11, 
November 2008. 

[32] Wright, M. W., White, T., and Mangini, N., “Data Par- 
allel Fine Relaxation (DPFR) Code User Manual Aca- 
dia Version 4.01.1,” NASA/TM-2009-215388, October 
2009. 

[33] Anderson, W. K. and Bonhaus, D. F., “An Implicit 
Upwind Algorithm for Computing Turbulent Flows on 
Unstructured Grids,” Computers and Fluids, Vol. 23, No. 
1, 1994, pp. 1-21. 

[34] Anderson, W. K., Rausch, R. D., and Bonhaus, D. F., 
“Implicit/Multigrid Algorithms for Incompressible Tur- 
bulent Flows on Unstructured Grids,” Journal of Compu- 
tational Physics, Vol. 128, No. 2, 1996, pp. 391-408. 

[35] Nichols, R. H., Buning, P. G., “User’s Manual for 
OVERFFOW 2.1,” University of Alabama and NASA 
Fangley Research Center, 2008. 

[36] Daso, E. O., Pritchett, V. E., Wang, T. S., Ota, D. K., 
Blankson, I. M., Auslender, A. H., “Dynamics of Shock 
Dispersion and Interactions in Supersonic Freestreams 
with Counterflowing Jets,” AIAA Journal Vol. 47, No. 
6, June 2009. 

[37] Codoni, J. R., Berry, S. A., “Analysis of Dynamic 
Data from Supersonic Retropropulsion Experiments in 
NASA Fangley’s Unitary Plan Wind Tunnel,” AIAA 
Paper 2012-2706, June 2012. 


Biography 


Guy Schauerhamer is an aerospace en- 
gineer for Jacobs Technology supporting 
NASA Johnson Space Center's Applied 
Aeroscience and CFD Branch. He re- 
ceived his B.S. and M.S. in Mechanical 
Engineering from Utah State University 
where his thesis was on six-degree-of- 
freedom moving mesh applications. His 
primary focus is in the field of com- 
putational fluid dynamics, in which he 
has been involved in the aerodynamic analysis of capsules, 
Ares V, full aircraft, the Morpheus lander, and Supersonic 
Retropropulsion. 


Kerry Zarchi (Trumble) is an aerospace 
engineer in the Aerothermo dynamics 
Branch at NASA Ames Research Center. 
Her primary focus is in the field of hy- 
personic computational fluid dynamics, 
in which she has been involved in the 
aerothermal analysis of Stardust, Orion 
Crew Exploration Vehicle, Space Shut- 
tle Damage Assessment, and Supersonic 
Retropropulsion. She received her B.S. 
in Mechanical Engineering from Michigan Technological 
University and her M.S. in Aerospace Engineering from Uni- 
versity of Minnesota, where her thesis was on flow control 
through laser energy deposition. 

Bil Kleb is an aerospace engineer in the 
Aerothermo dynamics Branch at NASA 
Langley Research Center. His main area 
of research has been the development 
and application of hypersonic compu- 
tational fluid dynamics codes such as 
LAURA and FUN 3D. He has a long 
history of providing aerothermo dynamic 
analysis of the shuttle and has con- 
tributed to numerous NASA space access 
studies such as X-33, X-34, X-37, and Orion. His current 
research focuses on verification and validation of compu- 
tational fluid dynamic models to enable prediction of Mars 
decelerator reentry performance. He has a B.S. and M.S. 
in Aeronautics and Astronautics from Purdue University, an 
M.B.A. from the College of William and Mary, and a Ph.D. in 
Aerospace Engineering from the University of Michigan. 

Karl Edquist is an aerospace engineer 
in the Atmospheric Flight and Entry 
Systems Branch at NASA Langley Re- 
search Center. His areas of expertise 
are computational aerothermo dynamics 
and aerodynamics, and he has authored 
or co-authored several papers in these 
fields. He is the Supersonic Retropropul- 
sion Lead for the Entry, Descent, and 
Landing Project as well as the Aerother- 
mal Lead for the Mars Science Laboratory entry capsule. 
He received his B.S. in Aerospace Engineering from the 
University of Colorado and M.S. in Aerospace Engineering 
from the University of Maryland. 






14 




